===== 描述性统计 =====
# A tibble: 9 × 7
  year  trait    mean    sd median   min   max
  <fct> <chr>   <dbl> <dbl>  <dbl> <dbl> <dbl>
1 2015  height 119.   65.4     112   3     390
2 2015  crow   150.   80.0     140   1     540
3 2015  diam    23.2  28.0      12   0.1   200
4 2020  height 115.   71.2     105   4     375
5 2020  crow   156.   91.1     145   5     528
6 2020  diam     7.72  9.27      5   0.2   145
7 2025  height 119.   65.6     110   3     380
8 2025  crow   166.   96.8     150   5     723
9 2025  diam    10.2  12.4       7   0.2   157

===== Kruskal-Wallis检验 =====
           Trait     chi_sq df       p_value
1         Height   7.483724  2  2.370991e-02
2 Crown diameter  10.893595  2  4.310085e-03
3 Basal diameter 463.259250  2 2.538236e-101

===== Dunn事后检验 =====
$Height
   Comparison          Z    P.unadj      P.adj
1 2015 - 2020  2.2587692 0.02389775 0.03584662
2 2015 - 2025 -0.1419989 0.88708090 0.88708090
3 2020 - 2025 -2.4943204 0.01261986 0.03785957

$Crown
   Comparison          Z     P.unadj       P.adj
1 2015 - 2020 -0.7156568 0.474203262 0.474203262
2 2015 - 2025 -3.1236408 0.001786284 0.005358851
3 2020 - 2025 -2.3768487 0.017461248 0.026191871

$Diameter
   Comparison         Z      P.unadj        P.adj
1 2015 - 2020 20.950996 1.837478e-97 5.512433e-97
2 2015 - 2025 15.115699 1.276100e-51 1.914150e-51
3 2020 - 2025 -6.710784 1.935814e-11 1.935814e-11


===== 异速生长模型系数 =====
# A tibble: 9 × 6
  year  intercept slope r_squared   p_value relationship
  <fct>     <dbl> <dbl>     <dbl>     <dbl> <chr>       
1 2015      0.574 0.832     0.739 0         logH~logC   
2 2020      0.215 0.891     0.774 0         logH~logC   
3 2025      0.966 0.738     0.544 1.08e-234 logH~logC   
4 2015      3.63  0.379     0.453 9.60e-153 logH~logD   
5 2020      3.69  0.541     0.622 2.13e-244 logH~logD   
6 2025      3.77  0.439     0.408 2.29e-157 logH~logD   
7 2015      3.68  0.452     0.605 8.35e-234 logC~logD   
8 2020      4.04  0.513     0.575 2.94e-215 logC~logD   
9 2025      4.16  0.405     0.348 7.94e-129 logC~logD   

===== ANCOVA交互作用检验 =====

--- Height ~ Crown * Year ---
Anova Table (Type III tests)

Response: logH
            Sum Sq   Df  F value    Pr(>F)    
(Intercept)   7.90    1   51.448 8.854e-13 ***
logC        395.80    1 2576.131 < 2.2e-16 ***
year          7.14    2   23.246 9.293e-11 ***
logC:year     7.33    2   23.860 5.068e-11 ***
Residuals   562.02 3658                       
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

--- Height ~ Diam * Year ---
Anova Table (Type III tests)

Response: logH
             Sum Sq   Df    F value    Pr(>F)    
(Intercept) 2807.13    1 11476.9996 < 2.2e-16 ***
logD         242.56    1   991.6995 < 2.2e-16 ***
year           2.67    2     5.4511  0.004327 ** 
logD:year     19.86    2    40.5932 < 2.2e-16 ***
Residuals    894.70 3658                         
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

--- Crown ~ Diam * Year ---
Anova Table (Type III tests)

Response: logC
             Sum Sq   Df   F value    Pr(>F)    
(Intercept) 2889.05    1 11910.588 < 2.2e-16 ***
logD         346.09    1  1426.811 < 2.2e-16 ***
year          29.60    2    61.010 < 2.2e-16 ***
logD:year      7.82    2    16.116 1.075e-07 ***
Residuals    887.29 3658                        
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
